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Abstract 

In an attempt to produce a reliable source of water, the water quality mod¬ 
el—QUAL2Kw—was put to through calibration and validation tests. This 
model was implemented in the river Tungabhadra of India. This worked 
quite well for the project at hand but showed certain irregularities. Different 
Quality Management methods exist currently and for this project, their im¬ 
pact in DO concentrations were examined based on: 1) The ability to aug¬ 
ment the flow; 2) The ability to oxygenate the source and; 3) The ability to 
modify pollution loads. With the impacts realized, the study helped us under¬ 
stand that local oxygenation is effective when considering the increment of 
Do levels. There was also the understanding that a combination of the fac¬ 
tors (flow augmentation, pollution load modification and source oxygena¬ 
tion) is a necessity as far as the minimum DO concentrations must be rea¬ 
lized. With this quality model results, there is confidence that the imple¬ 
mentation of the QUAL2Kw model will be an appropriate choice for future 
river water quality policies. 
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1. Introduction 

The run for clear river water qualities arises due to the increasing scarcity of us¬ 
able water in developing countries. In a major case, human factors have been 
found to be the primary factor that causes the unavailability of clear water hence 
the main objectives of the river basin management [1]. With the assessment and 
evaluation of human impacts on the quality of water as a major objective, there 
has been a major issue that breaches this objective. This issue is the prediction of 
the release of chemical substances into the river system by humans [2]. This im¬ 
pairs the efforts put in place to determine the impact of human activity in the 
aquatic ecosystem. 

Human activities generate massive amounts of contaminable materials on a 
daily basis. They could come from agricultural, municipal or industrial sources. 
The materials released could be large amounts of nutrients of organic materials 
and their introduction into rivers or streams accelerates eutrophication processes 
and reduces dissolved oxygen values. This is visible during low flow periods as 
the values are brought below threshold figures (Figures 1-7). The effect of this is 
an unbalanced aquatic ecosystem resulting in the death of aquatic life. This 
death could stem from extreme anaerobic conditions caused by low dissolved 
oxygen (DO) concentrations [3]. 


Study Area 
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Tungabhadra River Stretch 
Figure 1. Sampling sites along Tungabhadra river. 
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Figure 2. Mass balance equation in reach segment i. 



Figure 3. Model kinetics and mass transfer processes. The state variables are defined in 
QUAL2KW manual. Kinetic processes are dissolution (ds), hydrolysis (h), oxidation (ox), 
nitrification (n), denitrification (dn), photosynthesis (p), respiration (r), excretion (e), 
death (d), respiration/excretion (rx). Mass transfer processes are reaeration (re), settling 
(s), sediment oxygen demand (SOD), sediment exchange (se), and sediment inorganic 
carbon flux (cf). 
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Figure 4. Segmentation scheme to calculate K a river. 
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Figure 5. Calibration in Tungabhadra river for post-monsoon. 
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Figure 6. Validation in Tungabhadra river for pre-monsoon. 
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Figure 7. Dissoloved oxygen profile along river Tungabhadra with different BOD and 0.3 mg/1 TN limits at 
point load. 



In a general sense, any water quality management policy should strive to 
maintain pollution levels below threshold levels. This will ensure minimum DO 
concentrations on aquatic life. A minimum DO concentration is a vital aspect 
when considering aquatic life. According to the [4], the acute lethal limit of DO 
concentration for salmonids (type of fish) is about 3 mg/L. When considering 
the Tungabhadra river (which has varieties of aquatic life), the mild warm-water 
minimum can be put at about 4 mg/L—a day minimum. This is made to support 
less tolerant species and capricious fish inhabitants [5] [6]. The assimilative ca¬ 
pacity of a river should be kept at the optimal level in order to achieve the de¬ 
sired water quality [7] and it can be touched by: 1) The application of oxygena¬ 
tors [7], 2) The controlling of pollution loads [8] and 3) The controlling of river 
flow rates [9]. 

Strategies of water quality management are mostly a combination of “in¬ 
ter-disciplinary” complex resolutions. These decisions are based on assumptions 
to the changing controls of water quality [10]. The linkage between waste loads 
and resulting water qualities can best be shown via mathematical methods [11]. 
The most used mathematical model to evaluate the pollution of water bodies is 
the QUAL2E model [12]. This model was developed by the United States Envi¬ 
ronmental protection Agency (US EPA). While this has been the case, some li¬ 
mitations have been listed by users and examiners of the model [13] [14] [15] 
[16]. A major inadequacy of this model is the unavailability of parameters to 
handle aspect such as conversion of algal death to carbonaceous biochemical 
oxygen [17] [18]. 

Given the erroneous nature of the QUAL2E model, it was later modified. The 
modifications were inclusive of an expansion to accommodate computational 
structures, the addition of new constituent interactions such as de-nitrification, 
algal BOD and DO changes. These were caused by fixed plants. After this event 
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had occurred, Park and Lee in 2002 developed the QUAL2K model and in 2005, 
Pelletier and Chapra developed the QUAL2Kw model. The difference between 
this and the QUAL2K model (also developed by Chapra and Pelletier in 2003 
[19] was the intention to modernize the QUAL2E model. 

The QUAL2Kw model takes on new elements to optimize the water quality 
policy [20]. In the first instance, it makes use of two forms of carbonaceous bio¬ 
chemical oxygen to replace organic carbon. They are; slowly and rapidly oxidizing 
forms. Also, it aids the reduction of oxidation reactions to about zero at reduced 
oxygen levels by including anoxia and works explicitly when considering bottom 
algae attachments. The internal part of the model allows sediment-water changes 
of nutrients and dissolved oxygen, and its simulation is extended to accommo¬ 
date pH, quality of pore water (considering sedimentation) and de-nitrification. 

The QUAL2Kw model is a one-dimensional model that can be implemented 
on the Microsoft Windows platform. It is available on the internet and can be 
accessed via http://www/epa.gov/- . The QUAL2Kw is a steady state stream water 
quality model developed to simulate constituents including the pH, sediment 
oxygen demand, dissolved oxygen, ammonia nitrogen, organic phosphorus, total 
nitrogen, bottom algae, phytoplankton, total phosphorus, nitrite and nitrate ni¬ 
trogen, temperature, carbonaceous biochemical demand, ammonia nitrogen and 
organic nitrogen. 

The potentials of the model back the choice of the model as the framework for 
this study considering the Tungabhadra River. Within this study, I explain the 
application of the model and examine the impact of waste loads on water bodies. 
I also determine the utmost pollution loads that are acceptable by the stream 
without the violation of the least possible DO standards. 

2. Material and Methods 
2.1. Study Area 

River Tungabhadra is a part of the Krishna basin and formed by the combina¬ 
tion of the River Bhadra and River Tungha at Koodli which move towards down 
the eastern-incline of the Western-Ghats in the domain of Karnataka. The two 
waterways begin in MudigereTaluk of Chikmagalur District of Karnataka along¬ 
side the Nethrvathi, the Tunga and the Bhadra ascend at Gangamoola, in Vara- 
ha-Parvatha in the Western-Ghats at a height of 1198 meters (close Samse Vil¬ 
lage). From that point, the Tungabhadra walks through the fields to a distance of 
531 km. After intersection, Tungbhadra stream moves through Honnali, Hari- 
hara and Harapanahallitaluks of Davangere region. At that point, it moves 
through Hoovina-Hadagali, Hagari-Bomannahalli, Hospet and SiruguppaTaluks 
of Bellary area. In Siruguppa Taluk of Bellary location, it gets its tributary Veda- 
vathi River. The waterway shapes characteristic limit among Bellary and Koppal 
district and after that among Bellary and Raichur locale along its course. In the 
wake of entering Andhra Pradesh State, it moves through sanctuary town of 
Mantralaya and after that through Kurnool Taluk. Its rise from the ocean level is 
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539 m and the co-ordinates are 14.52”N and 75.8°E. Population of the Harihara 
city preliminary at 2011 was 85,600 (as per census) and the anticipated popula¬ 
tion for 2018 is 110,000. 

This study covers a 40.05 km stretch of river starting 18 km upstream of Ha¬ 
rihara town to about 22.5 km downstream. Stream water samples were Col¬ 
lected. from seven sampling-station viz. Mudenuru (SSI), Rajanahalli (SS2), 
Kumarapatnam bridge (SS3), Nalawagal (SS4), Nadhiharalli (SS5), Airani (SS6) 
and Heribedri and Kadathi (SS7), and Somalapura (SS8) with populations rang¬ 
ing from 45,000 to 55,000, are located on the bank of this river at the Upstream 
and downstream of Harihara town. For all these villages and Harihara town, this 
river is a primary source of drinking water. 

The river about two to three decades ago was in a drinkable state. It was 
deemed pure in both physical and ritual forms. This is however, no more feasi¬ 
ble. Ritual practices such as bathing, etc. are almost extinct due to the level of 
pollution in the river [5], The number of aquatic life in this holy river has also 
reduced drastically as it receives an enormous discharge of domestic and indus¬ 
trial effluent [21]. It is thus labeled biologically dead. It is because of occurrences 
such as this that made the selection of the 40.5km stretch of the Tungabhadra 
river feasible for this study. 

2.2. Data and Monitoring Sites 

Data for this study was collated between the 9 th and 19 th of June, 2018 as well as 
between the 2 nd and 3 rd of December, 2018. This counts as pre and post monsoon 
periods for the Tungabhadra river and its tributaries. The maximum number of 
hours spent was about 30. 

During the pre-monsoon season, sample events were set in order to monitor 
low flows at their critical state. Data pulled for this included air temperature, 
wind speed, relative humidity counted at meteorological data and they were ob¬ 
tained from the Department of Meteorology and Hydrology. Major points used 
for monitoring along the river include, sampling Station 1 (Ingalagondi and 
Mudenuru, chainage) at near to Jack well point, sampling Station 2 (Rajanahalli 
village) adjacent to Temple road, sampling Station 3 (Kumarapatanam bridge) at 
wastewater discharge from Harihara and Kumarapatanam town, sampling Sta¬ 
tion 4 (Nalavagalu) at near to HariharaPolyfiber discharge, sampling Station 5 
(Nadiharahalli) at just near from Nalavagal, sampling Station 6 (Airani) at just 
downstream of Temple, sampling Station 7 (Heribidri) at just downstream of 
confluence with sugarcane factory waste molasses and sampling Station 8 (So¬ 
malapura) at just downstream of temple. Observations were made between sta¬ 
tions Rajanahalli and Kumarapatanam at about 8.000 km chainage and the 
quantities of water abstractions noted were about 0.14 m 3 /s in pre-monsoon and 
0.36 m 3 /s in post-monsoon periods (Tables 1-5). 

2.3. Sampling and Analysis 

The parameters examined for water quality included; Temperature, pH, dissolved 
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Table 1. Water quality of surface water along Tungabahadra river in premonsoon. 


Sampling station 

Temp 

PH 

DO 

BOD 

TN 

TP 

Organic-N 

N0 3 -N 

NH„-N 

Organic-P 

Inorganic-p 




mg /1 

mg /1 

mg /1 

mg /1 

mg/l 

mg/l 

mg/l 

pg/i 

pg/i 

SSI 


7.2 

7.8 

2.44 

5.65 

0.27 

0.43 

0.19 

5.03 

184.1 

85.9 

SS2 


7.9 

8.2 

7.68 

11.48 

1.17 

2.465 

1.385 

7.63 

26.18 

905.8 

SS3 


7.65 

8.4 

5.62 

12.45 

1.15 

2.27 

1.88 

8.3 

876 

624 

SS4 

28 

8.2 

4.1 

7.6 

11.95 

2.15 

1.66 

1.89 

8.4 

1053.5 

1096.5 

SS5 


8.16 

4.85 

6.45 

14.23 

3.84 

2.42 

2.2 

9.6 

1167.4 

2627.5 

SS 6 


8.1 

5.5 

5.62 

15.32 

3.9 

3.62 

2 

9.7 

1095.9 

2804.1 

SS7 


7.8 

7.65 

4.4 

15.4 

2.85 

3.14 

1.96 

10.3 

883.5 

1966.5 

SS 8 


7.4 

8.85 

3.1 

13.2 

2.16 

2.44 

2.36 

8.4 

475.2 

1684.8 

Table 2. Water quality of surface water along Tungabahadra river in postmonsoon. 

Sampling station 

Temp 

PH 

DO 

BOD 

TN 

TP 

Organic-N 

NO 3 -N 

nh 4 -n 

Organic-P 

Inorganic-p 




mg /1 

mg /1 

mg /1 

mg /1 

mg/l 

mg/l 

mg/l 

pg/i 

pg/i 

SSI 


7.1 

7.6 

2.4 

5.05 

0.37 

0.6 

0.35 

0.95 

116.9 

253.1 

SS2 


8.1 

8.2 

7.8 

6.99 

2.6 

1.23 

0.61 

1.84 

1674.4 

925.6 

SS3 


8.43 

9.1 

7.5 

8.77 

3.48 

1.6 

0.74 

2.34 

1057.9 

1990.1 

SS4 

25.1 

8.12 

8.1 

6.3 

10.11 

4.48 

1.74 

0.68 

2.42 

1617.3 

2862.1 

SS5 


8.3 

6.16 

5.62 

11.84 

5.23 

2.65 

0.83 

3.48 

3059.55 

2170 

SS 6 


7.85 

6.34 

4.37 

9.8 

5.85 

2.09 

0.74 

2.83 

3591.9 

2258 

SS7 


7.7 

8.16 

3.2 

10.51 

6.1 

1.67 

0.45 

2.12 

1210.4 

4889.6 

SS 8 


7.4 

8.7 

2.6 

10.26 

4.15 

1.1 

0.76 

1.86 

938.315 

3211.68 


Table 3. River hydraulic characteristics for Tungabhadra river. 



season 

velocity 

depth 

flow velocity depth 


coefficient 

exponent coefficient exponent 

m3/s m/s 

m 

premonsoon 

0.408 

0.365 0.318 

0.459 

0.76 0.37 

0.28 

postmonsoon 

0.389 

0.373 0.112 

0.714 

5.51 0.73 

0.585 

Table 4. Composition (%) nitrogen and phosphorus at sampling sites along Tungabhadra 
river. 

sampling station 

organic-P 

Inorganic-P 

Organic-N 

NO r N 

nh 4 -n 

SSI 

49.88989 

50.11011 

9.745904 

5.146762 

85.10733 

SS2 

33.3188 

66.6812 

19.53435 

10.39561 

70.07005 

SS3 

53.28667 

46.71333 

18.23847 

11.76913 

69.59331 

SS4 

42.55022 

57.44978 

15.55095 

11.27096 

73.1781 

SS5 

44.45052 

55.54948 

19.69404 

11.23522 

69.24676 

SS 6 

44.75 

55.25 

22.47789 

10.30293 

69.41307 

SS7 

25.42131 

74.57869 

18.13962 

8.504455 

73.35593 

SS 8 

22.305 

77.695 

14.60305 

12.6431 

72.75385 

Average Values 

39.49655 

60.50345 

17.24803 

10.15852 

72.8398 
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Table 5. Calibration system parameters in the Tungabhadra river. 



Adopted 


Minimum 

Maximum 

Parameter 

values 

Units 

values 

values 

Carbon. 

40 

gc 

40 

40 

Nitrogen. 

7.2 

gN 

7.2 

7.2 

Phosphorus. 

1 

gP 

1 

1 

Detritus dry weight. 

100 

gD 

100 

100 

Chlorophyll-a. 

1 

gA 

1 

1 

Slow CBOD hydrolysis rate. 

3.83 

/day 

0 

5 

Slow CBOD oxidation rate. 

2.48 

/ day 

0.04 

4.2 

Fast CBOD oxidation rate. 

3.2 

/day 

0.02 

4.2 

Organic N hydrolysis. 

0.4 

/day 

0.02 

0.4 

Organic N settling velocity. 

0.1 

m/day 

0.001 

0.1 

Ammonium nitrification. 

9.53 

/day 

0 

10 

Nitrate de-nitrification. 

1.48 

/day 

0 

2 

Sediment de-nitrification 

0.94 

m/day 

0 

1 

transfer coefficient. 

Organic P hydrolysis. 

0.7 

/ day 

0.01 

0.7 

Organic P settling velocity. 

0.003 

m/day 

0.001 

0.1 

Inorganic P settling velocity. 

0.13 

m/day 

0 

2 

Sediment P oxygen attenuation 
half saturation constant. 

1.88 

mg0 2 /L 

0 

2 

Bottom Algae: 





Max growth rate. 

106.1 

mgA/m2/day or /day 

0 

500 

Respiration rate. 

0.31 

/day 

0.05 

0.5 

Excretion rate. 

0.43 

/day 

0 

0.5 

Death rate. 

0.008 

/day 

0 

0.5 

External nitrogen half sat constant. 

33.83 

AgN/L 

10 

300 

External phosphorus half sat constant. 

30.12 

AgP/L 

1 

50 

Light constant. 

40.76 

langleys/day 

1 

100 

Ammonia preference. 

71.81 

AgN/L 

1 

100 

Subsistence quota for nitrogen. 

6.06 

mgN/mgA 

0.0072 

7.2 

Subsistence quota for phosphorus. 

0.3 

mgP/mgA 

0.001 

1 

Maximum uptake rate for nitrogen. 

80.59 

mgN/mgA/day 

1 

500 

Maximum uptake rate for phosphorus. 

490.42 

mgP/mgA/day 

1 

500 

Internal nitrogen half sat ratio. 

2.06 

Dimensionless 

1.05 

5 

Internal phosphorus half sat ratio. 

1.88 

Dimensionless 

1.05 

5 


oxygen (DO), 5 days’ biochemical oxygen demand (BOD), Total nitrogen (TN), 
Total phosphorus (TP), Nitrite + Nitrate-nitrogen (NO s -N), Organic-nitrogen 
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(Organic-N), Ammonia-nitrogen (NH 4 -N), Organic Phosphorus (Organic-P), 
Inorganic Phosphorus (Inorganic-P), River flow (Q), River velocity (m) and 
River water depth (m). The float method was used to measure floatation. This 
was determined via the consideration of cross-sections and the measuring of ve¬ 
locities at 3 - 5 pints along river cross-sections. The parameters in this study 
were expressed in mg/L or pg/L. All parameters except pH, temperature (°C), 
velocity (m), depth (m) and flow (m 3 /s). Also, standard methods were applied 
during sampling preservation, transportation and analysis of water samples 
(APHA, 1989). During analysis for temperature and pH, the tools used were the 
portable thermometer and the pH-meter respectively. The remaining parameters 
were analyzed thus; DO was measured with a portable DO probe and in the la¬ 
boratory by titration. BOD and CBOD were determined by azide modification 
method (5 days’ incubation period at 20”C). TP was analyzed via digestion + 
ammonium molybdenum blue method. The samples in quote were digested with 
concentrated nitric + sulphuric acids before analysis of phosphate. This was ne¬ 
cessary to convert phosphates into orthophosphates form. TN was analyzed by 
the determination of Nitrite-N, Nitrate-N and Total Kjeldahl Nitrogen (TKN). 
Nitrite-N was analyzed by spectrophotometric method (diazo method). TKN 
was analysed using macro-kjeldahl method by following titration with standard 
hydrochloric acid. Ammonia nitrogen was analyzed by nesslerization method 
and by titration. Organic-N was analysed by the subtraction of the observed val¬ 
ue of ammonia nitrogen from TKN [22]. 


2.4. Modeling Tool 


With QUAL2Kw, there is the possibility to check the concentration ci within the 
water columns. This is because the model has a general mass balance equation. 
This though does not work for hyporheic concentrations. In using the general 
mass balance equation, parameters representing transport and loading are likely 
to be omitted while the column could be considered to be of reach /[ 20] [21]. 
During auto-calibration for the QUALK2w model, we get to employ the genetic 
algorithm (GA). This step is a necessity if we desire to achieve a well-placed 
goodness of fit for the model. The results obtained from this algorithm can then 
be compared with the measured data. When carrying out the comparison, we 
adjust the parameters used. These parameters are usually large in numbers [20], 
The determination of the best goodness of fit for the model can be done as a re¬ 
ciprocal for the weighted average of the RMSE. RMSE refers to the root mean 
squared error. It must be normalized before the weighted average can be consi¬ 
dered. 


^£l 
d t 


Q-i 


v. 


a 


■ab,i 


El 


El. 


W, 


c i +-jr( c i - 1 -c i )+-^(c M -c i )+ v 


+s : +- 


F’ 

-C'/v 


hyp,i 


(c 2 i -c,)(l) 


where Wj = the external loading of the constituent to reach i [g/d or mg/d], and 
Sj = sources and sinks of the constituent due to reactions and mass transfer me¬ 
chanisms [g/m 3 /d or mg/m 3 /d]. Exchange of mass between the surface water and 
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the hyporheic sediment zone is represented by the bulk hyporheic exchange flow 
in reach i [ E ' hyp; in m 3 /day] and the difference in concentration in the surface 
water (c) and in the hyporheic sediment zone (c 2 / ). 

2.5. River Discretization 

During this study, the total length of the Tungabhadra river (40.5 km) was dis¬ 
cretized into 81 reaches of 0.5 km each. The model then makes use of the head¬ 
water to define the boundary condition upstream. The segment of the river 
treated as headwater was the sampling station “1” while internal calculation was 
used to represent the boundary condition downstream. The hydraulic features 
for each of the sampling locations were determined using measured river geo¬ 
metries and velocities. With this, the model was able to allow the input of the 
river reach hydraulic characteristics (coefficients and exponents of velocity and 
depth) as empirical equations to estimate average water velocity ( V) and depth 
( D ) of the river. During the pre and post-monsoon seasons, flows, mean depths 
and velocities were measured and used to compute the coefficients a, c and ex¬ 
ponents b , d. 

Define abbreviations and acronyms the first time they are used in the text, 
even after they have been defined in the abstract. Abbreviations such as IEEE, SI, 
MKS, CGS, sc, dc, and rms do not have to be defined. Do not use abbreviations 
in the title or heads unless they are unavoidable. 

2.6. Input Data 

As earlier explained, the water quality input parameters measured included flow, 
temperature, pH, DO, BOD, organic nitrogen, ammonia nitrogen, nitrite + ni¬ 
trate nitrogen, organic phosphorus and inorganic phosphorus. Data on inorgan¬ 
ic suspended solids, conductivity, fast CBOD, phytoplankton, detritus and pa¬ 
thogen were however not measured. When carrying analysis, a default value was 
adopted. Other inputs to the model included wastewater, groundwater, river 
tributaries and abstraction. 

Majorly, the inputs required for the model come from the families of nitrogen 
and phosphorus. This set off data was not measured in the pre-monsoon season 
however the post-monsoon survey was considered in the extraction of these da¬ 
ta. According to [6], during base-flow in Melbourne (urban storm water), total 
nitrogen was found to consist Organic-N 36%, NH 4 -N 9%, N0 2 + N0 3 -N 39% 
and particulate Organic-N 16%. What this showed was that the nitrogen compo¬ 
sition did not vary between sites. The proof for this was in the amount of total 
dissolved nitrogen. It was found that it comes out to be Organic-N 42.9%, NH 4 -N 
10.7% and N0 2 + N0 3 -N 46.49%. 

A survey on international researches put the composition of total nitrogen at 
Organic-N 71% NH 4 -N 5% and N0 2 + N0 3 -N 24%. Organic-N was found to be 
a major constituent in international literatures (Organic-N) [6]. If we are to con¬ 
sider the Melbourne data in which particulate nitrogen accounts for 16%, Or- 
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ganic-N will account for 55%. The composition of total dissolved nitrogen is 
65.5% NH 4 -N 6% and N0 2 + N0 3 -N 28.5%. According to another survey per¬ 
formed by Kucuksezgin et al. (2005), the composition of total dissolved nitrogen 
and phosphorus are given to be Organic nitrogen 88.5%, NH 4 -N 3.7% and N0 2 
+ N0 3 -N 7.8%; Organic-P 81.1% and Inorganic-P 18.9%, respectively. In the 
Tungabhadra river however, the average compositions of TP were Organic-P 
39.5% and Inorganic-P 60.5% while the average compositions of TN were Or- 
ganic-N 17.24%, NH 4 -N 72.83% and Nitrite + Nitrate-N 10.15%. Given the dif¬ 
ferences in the compositions of nutrients, the proportions derived along the 
Tungabhadra river were used in the modeling. 

In this study, I used sets of equations to act as processes for the QUALK2w 
model. The processes in the model could be physical, biological or chemical. The 
parameters within these are usually numerous. Surveys were made in respect to 
the rating for usable parameters. Choices made were dependent on results ob¬ 
tained from sources such as the Environmental Protection Agency, the QUAL2Kw 
user manual (Pelletier and Chapra, 2005) and Documentation for the enhanced 
stream water quality model QUAL2E and QUAL2E-UNCAS [12]. A number of 
options could be considered when calculating re-aeration using QUAL2Kw 
model. For the purpose of this study, it is calculated as a function of river hy¬ 
draulics using the Owens-Gibbs formula [23]. The formula according to [24] can 
be used effectively when working with water bodies that have depths up to 0.4 - 
11 ft. They also are useful when working with water bodies that have velocities 
that range from 0.1-5 ft/s. the CBOD oxidation rate for the study was however 
placed within the range 0.04 - 4.2. 

2.7. System Parameter 

The physical, chemical, and biological processes stimulated by QUAL2Kw 
represented by a set of equations contain many parameters Table 5. The ranges 
of model rate parameters were obtained from various literatures including En¬ 
vironmental Protection Agency (EPA) guidance document [25], QUAL2Kw user 
manual [21] and Documentation for the enhanced stream water quality model 
QUAL2E and QUAL2E-UNCAS [12]. QUAL2Kw has eight options to calculate 
re-aeration rate as a function of the river hydraulics. We have used Owens-Gibbs 
formula [23], which was developed for streams exhibiting depths ranging from 
0.4 to 11 feet and velocities ranging from 0.1 to 5 feet/s [24]. The range of CBOD 
oxidation rate was assumed as 0.04 - 4.2 as in 36 rivers in USA [26]. The settling 
of CBOD is considered insignificant. The ranges of other parameters were as¬ 
sumed default as in QUAL2Kw. 

2.8. Model Implementation 

Data gotten from the post-monsoon period were used for calibration. The time 
step for this was left at 5.625 min. This aided the avoidance of instability within 
the model. During integration, Euler’s method was used and the goodness of fit 
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was done using different weights assigned to varying parameters. The weights 
considered were found to reduce the errors between measures and modeled pa¬ 
rameter values. The most influential parameter for the model was the weight for 
DO and it was given as 10. Other parameters such as TN, TP, temperature, BOD, 
pH were given values of 2 while others assumed 1. 

The model was looped until proper adjustment was obtained for the parame¬ 
ters. This process also aided the achievement of an agreement between the re¬ 
sults as well as the measurements. The population sixe used during the model 
run was 100 and the generation span was 50. According to [20] the population 
size of 100 offers better performance as that of 500. 

Without changing the calibrated parameters, the model was run on data sets 
different from the ones taken during the pre-monsoon period. This was neces¬ 
sary to test the ability of the calibrated model in order to effectively predict the 
water quality under varying conditions. After this had been concluded, the cali¬ 
brated model could then be used to experiment on water quality conditions 
during critical periods in the period. 

3. Result and Discussion 
3.1. Calibration and Validation 

The model calibration for the water quality was monitored from six locations. 
The results gotten from the calibrations are similar to measured data—with the 
presence of some irregularities. From these results, the standard error (between 
the observed and simulated values) are 11.2%, 18.3% and 13.55% for flow, depth 
and velocity respectively. The effectiveness and interrelationships between the 
values are expressed by R 2 and they are 0.95, 0.82 and 0.84 respectively. The rela¬ 
tive mean errors together with values for R 2 are given thus; 25.62% (0.979), 
31.1% (0.98), 44.4% (0.98), 30.5% (0.948), 24.13% (0.96) and 17%. these values 
represent DO, BOD, TN, TP, temperature and pH. After simulations, the results 
show that the water qualities of the Tungabhadra river do not meet the mini¬ 
mum DO standard (4 mg/L beyond 18 km). In the upper segment of the river, 
there is an indication of better water quality as the DO stands at 5 mg/L. This 
owes to the presence of sewers to manage domestic, industrial and rain waste- 
water. Collected wastewater can be treated and dumped back into the river at a 
different section—say 17.05 km. The dumped packages can be seen later at about 
13-17 km and 18-21 km. the positioning of the packages can be given as; 
Package 1: 13.5 km and Package 2: 7.05 km - 8.2 km. 

Past 12 km, the concentration of dissolved oxygen decreases gently till it 
reaches a minimum of 14.25 - 15.25 km. This is due to the presence of a high 
level of pollution (industrial waste). Re-aeration coefficient was found to be low 
due to the increase in the velocities. Also, past 25.25 km, there is a gentle in¬ 
crease in the concentration of dissolved oxygen. This is due to the increase in the 
addition of increased DO water from sampling station 5 and re-aeration. 

Aside these, there is also an increase in the concentration of total potassium, 
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total nitrogen and CBOD. This occurs 18 km due to discharge of local waste wa¬ 
ter drains and treatment plant effluent. The simulation shows that there is sig¬ 
nificant difference between the modeled and measured values. The difference is 
partly because of the non-inclusion of pollution effective materials such as hari- 
hara poly fiber and Grism industry. 

The validated results for the model have some exceptions. For example, the SE 
and R 2 values: 53.25% (0.98), 42.65% (0.93), 23.5% (0.978), 25.93% (0.95), 63.15% 
(0.95) and 22% for DO, BOD, TN, TP, temperature and pH respectively are be¬ 
tween the simulated and observed values. These values despite the differences 
are acceptable for developing countries when considering the financial resources 
for these countries and the higher accuracy of data analysis. 

3.2. Strategies for Water Quality Control 

Point sources also influence DO concentrations along the river. Point sources 
include wastewater treatment plants, local sewers and river tributaries. The in¬ 
fluence of the point sources was evaluated using a validated model during the 
pre-monsoon period. This was in difference to the 7-day consecutive low flow 
with a 10 year return frequency. During the evaluation, there was the assump¬ 
tion of the absence of water leakage from sewers and the consideration of the 
following: 

1) Pollution loads modification 

2) Flow Augmentation and 

3) Local oxygenation. 

3.3. Modification of Pollution Loads 

Throughout the duration of the study, point sources used were assigned fixed 
trial values. The values were CBOD as 15 mg/L, 10 mg/L, 5 mg/L and TN as 0.3 
mg/L. The point sources were from the Tungabhadra river and there was the as¬ 
sumption of modification for the waste water treatment plants. This was effected 
successfully due to the enforcement of policies in respect to this purpose. All 
dissolved oxygen did not however meet the minimum requirements of 4 mg/L. 

3.4. Augmentation of Flow 

In this research work we are assuming that augmenting the flow of pollution 
loads was feasible upon conclusion of water supply projects in Karnataka. For 
that project is estimated to supply about 3.1 m 3 /s of water to the Tungabhadra 
river. In addition, none of the dissolved oxygen profiles actually equal the re¬ 
quirements. 

3.5. Local Oxygenation 

The effects of oxygenators along critical points of the river produces strong oxy¬ 
genation through air entrainment. The amount of dissolved oxygen coming into 
a place can be figured out via the implementation of an empirical equation. The 
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equation relates the Do deficit above and below dam to the geometrical proper¬ 
ties of the weir, weir type, quality of water and water temperature. The weirs 
were used to evaluate the effects of the oxygenators. Upon the conclusions of 
simulations, three primary points were discovered. The points—12.5 km, 13.5 
km and 14.5 km—were found suitable for the installment of the weirs. The ap¬ 
proximate heights needed for the weirs were about 1.35 m at 12.5 km and each 
0.75 m at 13.5 km and 14.5 km locations. With this in place, it is now possible to 
meet the minimum DO concentrations. This can be done by limiting the flow of 
augmentation with 5 mg/L CBOD. This however does not work at point 12.25 
km in which DO concentration is 2.4 mg/L. The reductions in the concentration 
of dissolved oxygen observed at points 12.25 km, 13.25 km and 14.25 km can be 
sad to be due the effect of weirs at 12.5 km, 13.5 km and 14.5 km respectively. 
This by extension increases the depth of the water and inversely acts on the aera¬ 
tion coefficient behind the dams. 

4. Summary 

With data pulled during 2018, the QUAL2Kw model was calibrated and tested 
for the Tungabhadra River. The model was suitable for the task to be done al¬ 
though there were inconsistencies. During calibration, the SE and R 2 values be¬ 
tween the simulated and measured data for DO, BOD, TN, TP, temperature and 
pH were 25.62% (0.979), 31.1% (0.98), 44.4% (0.98), 30.5% (0.948), 24.13% 
(0.96) and 17%. Upon validation, these values were discovered to be 53.25% 
(0.98), 42.65% (0.93), 23.5% (0.978), 25.93% (0.95), 63.15% (0.95) and 22% and 
15% respectively. 

In this study, the QUAL2Kw model was implemented in order to simulate 
different strategies when considering water management. This was done based 
on some factors including; flow augmentation, local oxygenation and pollution 
loads modification. With the results obtained, it can be said that local oxygena¬ 
tion is effective in raising dissolved oxygen levels. Also, it is seen that a combina¬ 
tion of wastewater modification, flow augmentation and local oxygenation can 
be used to ensure the requirements for the minimum DO concentrations are 
met. This places the model as a top choice when thinking of river water policy 
options. 
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